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The self-consistent spin-wave theory is applied to investi- 
gate the magnetization distribution around the impurity in 
isotropic and easy-axis two-dimensional ferro- and antiferro- 
magents. The temperature dependences of host magnetiza- 
tion disturbance and impurity magnetization are calculated. 
The short-range order in the isotropic case is investigated. Im- 
portance of dynamic and kinematic interactions of spin waves 
is demonstrated. 



consistent spin- wave theory (SSWT). This theory was 
developed to describe thermodynamics of 2D systems 
and also successfully applied to quasi-2D jlQ 11 
and weakly anisotropic 2D magnets 1 1 1 . An important 
advantage of SSWT in comparison with the usual spin- 
wave theory is a qualitatively correct description of the 
strong short-range order above Tm- Besides that, intro- 
ducing slave fermions jl^] into SSWT allows to take into 
account kinematic interactions of spin waves and describe 
systems with not too low Tm values. In the following 
sections we treat the cases of different signs of exchange 
interactions in the host and between host and impurity. 
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I. INTRODUCTION 

In connection with extensive investigations of copper- 
oxide based superconductors, great attention is paid last 
time to studying magnetism of low-dimensional systems. 
Of particular interest is the problem of non-magnetic 
impurities in magnetic hosts. Numerous experimental 
results (see, e.g., Jl|-|j|) demonstrate that even small 
amount of substitution impurities (Zn, Fe etc.) in Cu02 
planes may influence strongly magnetic properties, e.g. 
lead to strong suppression of host magnetization. These 
facts have stimulated a number of theoretical works 
(see, e.g., In particular, the impurity problem 

for isotropic two-dimensional (2D) antiferromagnets at 
T = was investigated by the standard spin- wave theory 
||. However, the detailed consideration of the finite tem- 
perature situation is absent. Moreover, the usual spin- 
wave theory is obviously inapplicable, since this does not 
take into account adequately the short-range magnetic 
order which is a characteristic feature of low-dimensional 
magnets. 

On the other hand, the impurity problem for three- 
dimensional (3D) magnets was investigated within the 
standard spin- wave theory (see, e.g. 0). It was estab- 
lished that in the case of a weakly coupled magnetic 
impurity in a ferromagnet the standard spin-wave ap- 
proximation is insufficient already at T ~ Ti mp where 
Timp "C Tc is the energy of impurity-host coupling. 
Inclusion of dynamic and kinematic interaction of spin 
waves within the Tyablikov approximation j(| leads in 
this case to occurrence of an anomalous temperature de- 
pendence of impurity magnetization. Therefore it is in- 
teresting to investigate the impurity problem for two- 
dimensional (2D) systems, such as ferro- and antiferro- 
magnets (FM and AFM) with small anisotropy or in- 
terlayer coupling, which are required to produce a finite 
value of the magnetic ordering temperature Tm ■ 

In the present paper we consider weakly anisotropic 
2D impurity magnetic crystals with the use of the self- 



II. FERROMAGNETIC IMPURITY IN 
FERROMAGNETIC HOST 

The Heisenberg Hamiltonian of a FM crystal with a 
quadratic lattice, containing a ferromagnetically coupled 
impurity at the site i — 0, reads 



where 



H A = -D^SD 2 - W^SfS. 



(1) 
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is the Hamiltonian of the easy-axis anisotropy, H is ex- 
ternal magnetic field. In the nearest-neighbor approxi- 
mation the non-zero exchange integrals are 



J', j = 0ori + £ = 
J, i, i + S ^ 



(2) 



where S denotes nearest neighbors, J > 0, J' > 0. 

Following to Ref. [O we use in the FM case for i =/= 
the representation |12| 



5+ = V2Sa, , Sf = S- a\a t - (25 + l)c\ Cl 



(3) 



,t„t. 



S i = V2S(a\ - ) 



2(25 + 1; 



25 



where a\ , ai are the Bose ideal magnon operators and 
cj , Ci are the auxiliary Fermi operators at the site i which 
take into account the kinematic interaction of spin waves. 
For i = one has to replace in (||) 5 — > 5' with 5' the 
impurity spin. Note that in the paper || only the Bose 
operators were introduced. 

To satisfy the condition Si(H = 0) = in the param- 
agnetic phase we introduce the Lagrange multipliers ^ at 
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each lattice site, which corresponds to the constraint of 
the magnon occupation number at T > Tq. These multi- 
pliers play the role of a local "chemical potential" for the 
boson-fermion systems. Introducing fa permits to correct 
the drawback of the standard spin-wave theory which is 
inapplicable at T > Tc since the magnetization formally 
becomes negative. Unlike the approach of Ref . || , we do 
not assume ad hoc the condition S — in the ordered 
phase where we have p,i = 0. Thus in our approach the 
magnon number is not conserved at T < Tq and the Bose 
condensation || does not take place. However, it may be 
shown that the results of both approaches are identical 
at low temperatures where kinematical interactions are 
small. 

Further we perform decouplings of the quartic forms 
which occur after substituting (g) into (|l|). Introducing 
the averages 



£i,t+<5 — S{ + s+ < a-idi+S > 



(4) 



we derive the quadratic Hamiltonian of the mean-field 
approximation 

H = X! &,i+sJid+s ajai - a\ +s a t + (25 + l)c]cj (5) 



\a\a i + {2S+l)c\c i 



Ha 



This Hamiltonian differs from that of the standard spin- 
wave theory in two points. First, the averages are 
introduced which take into account the dynamical inter- 
action of spin waves in the lowest Born approximation 
(see below). Second, the Fermi operators enter to ac- 
count the kinematic interactions of the spin- waves. 

Following to Ref. |l3|, we treat the influence of the 
magnetic anisotropy by neglecting quartic forms in Ha 
to obtain 

h a = -h a J2s? 

i 

= -H A J2[ S - a l a * - ( 2S + l ) c \°i 



with the anisotropy field Ha 



H a = (2S - 1)D + 



(6) 



Note that effects of the true magnetic field H and the 
field Ha are different in the paramagnetic phase since 
chemical potentials fii are calculated at H = 0, not Ha = 
0. Thus the field H yields a finite magnetization at any 
temperatures, whereas the field Ha induces a slight shift 
of Tq only. Thus the field Ha describes correctly the 
effect of the easy-axis anisotropy. In the limit Ha <C 
J under consideration effects of single-site and two-site 
anisotropy are the same, although concrete expressions 
for the field Ha in (m are different. 



For an ideal crystal ^^+5,^ do not depend on i and 
the diagonalization of the Hamiltonian (|^) is easily per- 
formed §,0. At the same time, for the impurity system 
this is a complicated task since the unknown dependence 
S,i,i+8, Hi which is to be determined self-consistcntly. 
However, as follows from the below calculations, Ci,i+s 
and practically coincide with the corresponding quan- 
tities for the host, £m and /z, except for nearest neighbors 
of impurity. Also in the FM phase as a function 

of i varies slower than the magnetization, and Hi = 0. 
Therefore we may put in 



'-.i.i-< > ■ = <1 i + S = , m - h 

£m i otherwise 



Sho, i = 
5 fix, i + S = 
0, otherwise 

(7) 



Note that £ ^ £' because of non-Hermiticity of the rep- 
resentation (p|) . Taking into account (|^) the spin correla- 
tion function of impurity spin with its nearest neighbors 
has the form 



K 



(8) 



Under the approximation [?] the Hamiltonian (^|) takes 
the form 



H = Ho 



V 



where 



Ha = J£m [4 a i - 4+s a * + ( 2S + l ) c U 
■18 

+ {H A +H-fx)J2 [a\ai + (2S + 1) 



(9) 



(10) 



is the standard SSWT Hamiltonian without impurities 
10 and 

V = (J'£ - J£ M ) [°o°o - a l a ° + i 2S ' + !)4co] (11) 

<5 

+ (j'c - j£ M ) [ a l a s - a l a s + ( 2S + l)4 c a 

+Sfi blb + 5/ii^oJo 4 
s 

To diagonalize H we introduce the Green's functions 

G°» = « a>j »°= £ -I—e***-*) (12) 

q H 

Gij(w) = < aj\a\ > w 



where the index means that statististical averages are 
calculated with H , 

Eq = ^m(Jo-Jh)+Ha+H- h, Jq = 2 J(cos q x + cos q y ) 



2 



In the limit i? > 1 we find by using the saddle point In the case of the pure system (V — 0) we have G — G° 
approximation (see, e.g., §) and the values Si, £,i.i+s, Mi are independent of i, so that 

the system of equations (p7|), ( Eg ) reduces to 
„ i A J exp(zvW^"^)/^ 1/4 ^ 1/2 1 « ("/J) 1/2 R, u> « 1 
Gor[U}) ~ \ - ln(w/ J) (a;/ J) 1 /^ « 1 

(13) * * " ' J ° 

The perturbation V can be written in the matrix form 



£ = = 5 + ^Vj k iV(£; k ) (19) 

3 = 5-^ iV(S k ) + (25 + 1)JV(£?/) 



i„3=0 



(14) 



i=0 



where the indices i, j enumerate the impurity site and its 
four nearest neighbors. From (Jll|) we have 



V = 



/ 4e 7 7 7 7 \ 

i p ' 

i p 

7' p 

V 7' p J 



, R=(2S + 1) 



P 
P 
P 



(15) 



7' = J'e- £ = 7' + <W4, 

7 = J'€'-J£m, p = i + 5pi 

Then we have the expression for the perturbed Green's 
function [pi: 



G{u) = [1-G°(w)V]- 1 G°(uj) 



(16) 



where G(uj), G°(ui) are submatrices of matrices 
Gij(uj), G®j(uj) with i, j — 0...4. Further we calculate the 

matrix G from dl6| ) and the averages < a\a,j > from the 
spectral representation. Then we derive from (||) , (0) the 
system of self-consistency equations 



+00 

£ = ^1 + / ^iV(w)ImG 10 H, 



(17) 



e = s a 



dto 



N(u})ImG 01 {u/) 



where N(ui) = l/(exp(u>/T) — 1) is the Bose distribution 
function. The integration region in ( |i"7j ) is in fact a < 
w < ^mJo + a, a = Ha + H — p,. The expressions for 
the site magnetizations take the form 



+ OC 



S = S'- I —N(Lu)lmG o(u) + (2S' + l)N{E ) 



— OO 
+ OO 



Si = S- I —N{u)TmG u (u) + (2S+l)N(E 1 ) (18) 

7T 



where Ef = (2S + 1) J^m + ol. One can see that £ de- 
pends on temperature due to dynamic magnon-magnon 
interactions; at low temperatures the corrections are pro- 
portional to T 5 / 2 , as well as in the Dyson's theory 

The results of numerical solution of the equations Q19|) 
at different values of Ha are shown on Fig.l. The results 
of numerical solution of Eqs. ( |l9[) for different Ha are 
shown on Fig.l. While the magnetization is strongly de- 
pendent from the value of Ha, the dependence £(T) is the 
same to calculation accuracy at an arbitrary Ha/ J <C 1. 
This dependence coincide with those for the ferromag- 
net in earliar variants of SSWT || at low temperatures. 
However, at T ~ J the short-range order parameter £ 
demonstrates a sharp decrease rather than vanishing. 
Thus introducing the Fermi operators removes the un- 
physical transition with vanishing of short-range order 
parameter. At finite values of Ha the value of the Curie 
temperature is finite. At small Ha <C J we have (cf. [fflj) 



Tn = 



4irJS 2 /ln(T/H A ) 
4nJS 2 /ln(ir 2 J/H A ) 



1 < \n(J/H A ) < 2irS 
ln{J/H A ) > 2-kS 



(20) 



Now we turn to the consideration of the impurity sys- 
tem. The results of numerical calculations of magnetiza- 
tions So, S\ vs. temperature according to ([l7|) , (|is|) in the 
zero magnetic field are presented in Figs 2,3. In Fig. 2 the 
results of the standard spin- wave theory (SW) which cor- 
respond, in our notations, to £m =£ = £' = S, fermion 
occupation numbers N(Ei) being replaced by zero, and 
the spin-wave theory with introducing fermions (SWF) 
are also presented for comparison. We see that the im- 
purity magnetization has an anomalous behavior at tem- 
peratures T ~ J' '. On the inset on Fig. 2 this dependence 
is shown at a different Ha/ J, J' /J- The sharp decrease 
of impurity-site magnetization at T ~ J' can be eas- 
ily obtained already in the simple mean-field approxima- 
tion, however the detailed description of the this behav- 
ior requires a more complicated methods. The standard 
spin-wave, as well as the SWF solution do not show this 
anomaly, so we can conclude that it is caused by both 
dynamic and kinematic interactions of spin waves. The 
situation is similar to the 3D case where using the Tyab- 
likov approximation results in a strong modification of 
the magnetization behavior in this temperature interval 



where E l = (25,; 
at the site i. 



1) J£m + a — 8 pi is the fermion energy 



In the ground state the disturbance of magnetization 
is localized at the impurity site and equals to S' — S. To 
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calculate the magnetization distribution around impurity 
at finite temperatures we need the full matrix G. It may 
be shown (see, e.g., that the latter quantity is given 
by 



G = G U + G 0N V 



1 



1-G°V 



-G 



No 



(21) 



where G 0N is the submatrix of with i = 0..4, j = 

Q..N, and G N0 is the conjugated matrix. Using (|2l| ) we 
can find the averages needed. The results of numerical 
calculation of magnetization disturbance for different val- 
ues of J'/ J, Ha/ J are presented in Fig. 4. One can see 
that at R > all results are practically the same, this 
takes place also in the limiting case with J' = (or in 
the case of vacancy with S' = 0). One can see that the 
change of magnetization around the impurity rapidly de- 
creases with increasing distance from the impurity site, so 
that the magnetization disturbance practically vanishes 
at the distance of 4 coordination spheres. 

In the 2D isotropic magnets where the long-range or- 
der at finite temperatures is absent we have to treat the 
short-range order parameters £y only. The chemical- 
potential corrections Sfii, i = 0, 1 are defined from the 
condition 



S 



S- < a\ ai > -(2S + 1) < c\c, > 

+ oo 



H=0 



(22) 



JV(w)ImG ii (w) 



H=0 



+ (2S+l)N(E t )\ H=Q =0 

In zero magnetic field the solution to eqs (|l7|), ( ^8|) 
is not unique. The absense of the unique solution in 
the paramagnetic phase is apparently the shortcoming of 
our approach. Since at small values of H the solution 
is unique, it is natural to take H to be small, but finite 
(in numerical calculations we have taken H = 0.005J). 
Although the site magnetizations are changed strongly 
with changing H (the susceptibility x = dS/dH is diver- 
gent near T = 0), it may be checked analytically that the 
derivative d^/dH remains finite at H — > 0, so that the 
values of £, £' weakly depend of H. This may be verified 
also by numerical calculations. 

The numerical procedure is as follows. To find the 
short-range order parameters £, £' we solve the system 
of equations ([jj]). At each iteration for given values of 
the corrections to the chemical potential 5(j,Q,5fj,i 
and mag netizations So, S\ are determined from eqs J22] ) 
and (jig), respectively. Results of numerical calculation 
of the short-range order parameters and correlation 
function K for the case of a weakly coupled impurity for 
the 2D ferromagnet are shown in Fig. 5. One can see that, 
owing to sharp decrease of the correlations between 
the impurity site and its nearest neighbors decrease with 
temperature more rapidly than those in an ideal crystal. 



III. ANTIFERROMAGNETIC IMPURITY IN 
FERROMAGNETIC HOST 

Further we consider an AFM impurity in FM host 
(J > 0, J' < in ||). After passing to the local coor- 
dinate system at the impurity site, we have to use the 
representation 

S+ = V^bl , = -S' + blb + (2S' + l)4d (23) 

1 2 S" 4- 1 

So = V2V(b - —blb b ) 2-j±-4d b 

where &J, b are the Bose operators, dj, do are the Fermi 
operators. Then, in the mean-field approximation, the 
Hamiltonian (|l|) takes the form 

W = ij ^2 £m+<5 a\a,i - a\ +s ai + (2S + l)b\bi 

i, i+S^O 

+i Jl i J2 [ a l as - b ° as + t 2S ' + v b l bs 

s 

H' [bX - b\a\ + (2S' + l)cjc ] } (24) 
+ ^{H A +H- IH ) \a\a t + (2S + l)b\b, 

+ (H A + H - mo) [blbo + {2S' + l)c c 
where 

£ = S + < 4al > 
£' = Si+ < a s d > 

As in ferromagnetic case, we use the approximation 
£i,i+S — S.M + 6 £ 0). To diagonalize (24) we in- 
troduce, following to H, the "hole" creation and annihi- 
lation operators aj, ao by the canonical transformation 



a 



-d 



As well as in the case of FM impurity, we use the approxi- 
mation (Q). We introduce also the Green's functions jl2] ) 
and represent the Hamiltonian as (^) with the parameters 
of the matrix V ( |l5| ) 



-- -J'C-J^M+SfJ.o/4,, 
J'i - J£ M + 6m 



Then we have the same equation ( |16| ) for the full 
Green's function as in the case of FM impurity, the self- 
consistency equations also has the same form (^7|), (|Tg|). 
Unlike the FM impurity case, the full Green's function 
has a pole at to — —ujo < ||. To take into account the 
contribution from this pole to the averages needed we de- 
form the integration path in the spectral representation 
for the Green's function in the complex plane: 



4 



< a] a, >= - / dujN(u;)Im.Gij(u>) (25) 

— OO 

= / ^7N(u)Gij(uj) - TGij(0) 



n 



The contour C is selected in such a way that all sin- 
gularities of G(u>) lie inside C, but all the frequencies 
uj n = 2imT (n ^ 0) lie outside it. The last term in ( p^ ) 
corresponds to the contribution from uj = which is to 
be subtracted explicitly. 

Results of numerical solution of Eqs. (jl?]), ( |l8| ) with 
using of fl25| ) for different values of impurity-host coupling 
and H = are presented in Figs. 6, 7. 

The AFM impurity induces the disturbance of host 
magnetization already at T = 0. Using the sum rule 

7T Vim + itf) - G&(w + W)] = ^-Imlndct[l - G (uj)V] 



5+ = V2Sb\ , Sf = -S + b\b l + (2S + l)4<k (29) 
Sr = V2~S( bi - ±b\bA) ~ 2 -^p±d\d i b i 

for i G B where a\,a i and 6^,^ are the Bose operators, 
c\,Ci and d\,di are the Fermi operators. After standard 
decoupling the Hamiltonian takes the form 



a l a i 



\+S a i 



{2S 



b l+sh+s - b]+s a t 



\Ji,i+8\ £,i,i+8 

+ Y,(Ha + H- [ojoj + (25 + l)c|c 



(30) 



+ (25+l) C J Cl 



^ + (25+1)^ 



where 



(26) 



ffA = (25-l)D + s5^|»h ii+tf | 



(31) 



which follows from ( pl| ) and taking into account that 
dct[l — Go(u)V] has a zero at uj = — ujq we obtain 



< b£bo >t=o= E < 



afcn >t=o 



(27) 



and 



<b ] - f 



J i+S u 'i 



a] >, £i,i+s = Si+ < a ib , 



i+5 



> 



i>0 



(32) 



so that the total disturbance of magnetization equals to 
S + S' |||. The distribution of magnetization around the 
impurity site is shown in Fig. 8. At large R the contri- 
bution from the pole uj — —ljq gives main contribution 
to the magnetization disturbance which is proportional 

to exp \ —R-\/uj~q/~J\ /R and differs from that in the 3D 

case B by preexponential factor only. 



IV. THE CASE OF AN ANTIFERROMAGNETIC 
HOST 

Now we consider an antiferromagnet with the Hamil- 
tonian 



For the correlation function K we have the same expres- 
sion (||) as in the FM case with £ = £oi>£' = Coi- To 
diagonalize PC]) we introduce the operators 



_ J a,i i e A 
i€B 

and the Green's functions: 
G«H =<^ AAA] ^> u = GlAQ) x 



r. 



i,j e A 
i,j € B 
otherwise 



(33) 



n 



where 



with Jij < 0, rjij < 0, Q =(n,ir), Hq is the staggered 
magnetic field. In the case of two sublattices A, B and 
for the antiferromagnetically coupled impurity spin in the 
A sublattice we have to use the representation 



r = 
A = 



X + UJ 



1/2 



A — uj 

Jo |€m + h a + h 



, CI = A - \/A 2 - uj 2 , 



S+ = V2Sa t , S* = S~ a\ ai - (2S + l)cjq 

rr— n&t t iff, 2(25+1) t 
S i =V2S{a\- —alalaA — a\ 



(28) 



for i S A and 



Using the approximation (Q) we get the same expres- 
sion for the Green's function ( |l6| ) with G°(w) being the 
5x5 submatrix of G^(uj), and analogously for G(uj). In 
this designations the self-consistent equations for the site 
magnetizations and short-range order parameters has the 
same forms as in FM case (|l8|), (|l7]). In the case of the 
pure system we now have G = Go so that 
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2 Jo ^ ££F 



coth 



AF 



2T 



s = s-J2 



^ Jk coth Ek 



AF 



k 2£jf 



2T 



+ (25 + 1)^(25,) 



where 



^k F = V(\Jo\dM + a) 2 ~(MM) 2 



The results of numerical calculations for the AFM im- 
purity system case are shown and compared with those 
for the FM case in Fig. 9. In the case of an impurity 
spin, which is weakly coupled to the host, the behavior of 
magnetization in AFM and FM situations is very close, 
except for the region near the magnetic ordering tem- 
perature (TV > Tc because of quantum fluctuations). 
At the same time, the nearest-neighbor magnetizations 
are strongly different and demonstrate a behavior, typ- 
ical for the corresponding hosts. One can also see that 
SSWT leads to unambiguous result at low T, where the 
impurity magnetization turns out to be greater than the 
host one. The difference between magnetizations of im- 
purity and host increases with lowering value of J'/ J and 
decreases with increasing temperature. Thus SSWT pre- 
dicts strong influence of quantum fluctuations on magne- 
tization in the case of weakly magnetic impurities. 

The results of calculating the short-range order param- 
eters £, £' and the correlation function K in the isotropic 
case Ha = are presented in Fig. 10. We use the same 
procedure as in the FM case. The parameter £ has a non- 
monotonic temperature dependence. At the same time, 
the temperature dependence of the correlation function 
of the impurity spin with its nearest neighbors is mono- 
tonic and more rapid than that for correlation functions 
between spins in the host. 

To calculate the total magnetization disturbance we 
use the sum rule for the Green's functions (B31) 



^V(-l)Tm Gu(uj + iS)-G%(co+iS) 



d_ 



Imln det[l 



(34) 



Since det(l — G°V) has no zeros at oj < 0, we obtain at 
T = 



SM = S'-S - -Im Indetfl - G°V) 

7T 



s' -s 

(35) 



This result is valid also for a vacancy if we put 5" = 0. For 
a ferromagnetically coupled impurity we have to replace 
in © S' -» -S'. 

The distribution of magnetization around impurity in 
the ground state of a 2D isotropic antiferromagnet is 
shown in Fig.ll. The magnetization of each sublattice 
decreases, so that corrections to the host magnetization 



have alternating signs. The values of sublattice magne- 
tization disturbance are close to those in the spin-wave 
theory ||. At large R, main contribution to the distur- 
bance of sublattice magnetization comes from frequencies 
lu -c J. Expanding (|2l]), ( ]33| ) up to first order in uj/J we 
derive 



S < A^At >~ l/Ri 



(36) 



Note that in the 3D case this quantity demonstrates a 
more rapid decrease 1/i? 4 , which may be obtained in the 
same manner. 



V. CONCLUSIONS 

To conclude, we have investigated 2D magnets with 
impurities for different signs of exchange integrals within 
the framework of self-consistent spin- wave theory 
This theory permits to calculate both magnetization dis- 
tribution and the correlation functions (short-range order 
parameters). For T — modifications of the results of 
the standard spin-wave theory are small. At the same 
time, for finite temperatures, corrections owing to dy- 
namic and kinematic interactions of spin waves turn out 
to be important. It should be stressed that despite the 
absence of long-range order in the isotropic 2D magnets 
at T > 0, the temperature dependence of the imurity- 
host correlation function K (g) is similar to that in the 
3D case, although in the latter case main contribution to 
K equals to SqSs- 

The distribution of magnetization in the ground state 
was investigated in detail. In the nearest-neighbor ap- 
proximation considered, the host magnetization distur- 
bance decreases rapidly with distance from impurity, and 
the total change of magnetic moment equals to — S ± 5" 
depending on the sign of J' . More interesting situations 
occur in the case of the long-range exchange. So, in the 
case of FM impurity in the FM host with sufficiently 
- :<sjftjin^[p5|2gative next-nearest impurity- host exchange J", 
the total magnetization change equals 



SM = S'-S- 2z 2 S 



(37) 



with Z2 the corresponding coordination number. In the 
case of FM impurity in the AFM host with large positive 
J" we have 



SM =S'-S+ 2z 2 S. 



(38) 



It would be of interest also to investigate the problem 
of a current carrier in the AFM host within a similar 
approach (e.g., within the t — J model, cf. 0). 

The work was supported in part by Grant No. 95-056 
from the State Scientific and Technical Program "Actual 
Researches in Condensed Matter Physics", subprogram 
"Superconductivity" . 
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Figure captions. 

Fig.l The temperature dependences of the magne- 
tization for the pure 2D ferromagnet with S = 1/2, 
H A /J = 1CT 2 (solid line), H A /J = 10~ 3 (dashed line) 
and short-range order parameter £ (dot-dashed line which 
is the same for all three cases: Ha/ J = 0, 10~ 3 , 10 -2 ). 

Fig. 2 The temperature dependence of the magneti- 
zations for impurity site So, and for nearest-neighbor 
sites S 1 (S = S' = 1/2, H A /J = 1(T 3 , J'/J = 0.15). 
The results of the magnetization for the impurity site 
in the standard (non-self-consistent) spin-wave approach 
without (SW) and with (SWF) introducing Fermi op- 
erators (see (||)) are presented for comparison. On the 
inset the temperature dependence of the impurity site 
magnetization So at Ha/ J = 1CP 3 , J'/J = 0.15 (solid 
line), H A / J = 10~ 3 , J'/J = 0.05 (short-dashed line), 
H A / J = 10~ 2 , J'/J = 0.15 (long-dashed line). 

Fig. 3 The temperature dependence of the short-range 
order parameters (solid lines, left scale) and corre- 
lation function K — (S0S5) (dashed line, right scale) for 
the same parameter values as in Fig. 2. Arrow shows the 
value of the Curie temperature. 

Fig. 4. The distribution of magnetization around impu- 
rity for the same parameter values as in Fig. 2, T = 0.3 J. 
Arrows show the value of magnetization disturbance at 
the impurity site (R = 0). 

Fig. 5. The temperature dependence of the parameters 
and spin correlation function K = (SqSs) in the 
isotropic 2D ferromagnet with S = S' = 1/2, J'/J = 
0.15. For comparison, the corresponding short-range or- 
der parameter in the ideal crystal, £m, is shown. 

Fig. 6. The temperature dependence of the magnetiza- 
tions for impurity site So, and for nearest-neighbor sites 
Si in the case of antiferromagnetic impurity in the ferro- 
magnetic host with S = S' = 1/2, H A / J = 10~ 3 , J'/J = 
-0.15. 

Fig. 7. The temperature dependence of the magnetiza- 
tions for impurity site So, and for nearest-neighbor sites 
Si in the case of antiferromagnetic impurity in the ferro- 
magnetic host with S = S' = 1/2, H A / J = 10~ 3 , J'/J 

I 

Fig. 8. The distribution of magnetization around anti- 
ferromagnetic impurity in the 2D isotropic ferromagnet 
at T = 0, J'/J = 0.15 (solid line), J'/J = 0.05 (dashed 
line). 

Fig. 9. The temperature dependence of magnetizations 
for impurity site So, and for nearest-neighbor sites Si, 
in an antiferromagnet with S = S' = 1/2, Ha/ J = 
10~ 3 , J'/J = 0.15. Dashed lines show the correspond- 
ing results for a ferromagnet (Fig. 2). 

Fig. 10. The temperature dependence of the short- 
range order parameters £, £' and correlation function 
K = — (S0S5) for the same parameter values as in Fig. 9. 
Arrow shows the value of the Neel temperature. 

Fig.ll. The distribution of magnetization around im- 
purity in the 2D isotropic antiferromagnet at T — 0. The 
inset shows the picture at large R. 
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